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ABSTRACT 

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission- 
line galaxies at z '--^ 2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with > la 
detections of emission lines to / > 2.5 x 10~^® erg s~^ cm~^, means that the galaxies in the sample 
are typically ~7 times less massive (median = than previously studied z ^ 2 emission- 

line galaxies. Despite their lower mass, the galaxies have [Oiii]/II/3 ratios which are very similar to 
previously studied z ~ 2 galaxies and much higher than the typical emission-line ratios of local galax- 
ies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the 
two-dimensional spectra of the galaxies for this purpose. In the stacked data the [O ill] emission line is 
more spatially concentrated than the H/3 emission line with 98.1% confidence. We additionally stack 
the X-ray data (all sources are individually undetected), and find that the average L[oiii]/-^o.5-iokoV 
ratio is intermediate between typical z ^ obscured active galaxies and star-forming galaxies. To- 
gether the compactness of the stacked [Oiii] spatial profile and the stacked X-ray data suggest that 
at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei. 



Subject headings: galaxies: abundances 
evolution 



galaxies: active — quasars: emission lines — galaxies: 
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1. INTRODUCTION 
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Rest-frame optical spectra contain a wealth of informa- 
tion about galaxy properties. In particular, the strengths 
and ratios between coUisionally excited and recombi- 
nation emission line s can be used to infer star forma- 
tion rate (SFR, e.g. lKennicuttlll998[ ). gas-phase metal- 
hcity (e.g.lKewlev et al.l 20011) and nuclear activit y (e.g. 
[Baldwin. Phillips fc Terlevichlri98lt IKewlev et an i2006) . 
However applying these techniques to galaxies at higher 
redshifts is challenging for two reasons: (1) at z > 1.5 
optical emission lines are observed in the near-infrared 
(NIR), where ground-based spectroscopy is challenging, 
and (2) the characteristics of z > 1 galaxies may be phys- 
ically different from local galaxies. 

The Wide Field Camera 3 (WFC3) G141 slitless 
grism on the Hubble Space Telescope (HST) addresses 
the first problem, with far improved capabilities for 
NIR spectroscopy tha n ground-based observatories (e.g., 
iStraughn et"aLl 120111) . Most of the previous ground- 
based studies of z > 1.5 galaxies hav e been biased 
to massive and UV-luminous systems fShaplev et al 
2004; Erb et al. 2006; Hayashi et al. 2009; Onodcra et a! 
|2010( ). probably missing most of the more numerous, less 
massive galaxies at 1.5 < z < 3. Only a handful of 
lower-mass galaxies at z '-^ 2 have been studied, revealing 
much lower metallicities than both z ~ galaxies of the 
same mass and th eir massive and eyolved counterpa.rts at 
the same redshift (lYuan fc Kewlevl[2009l : lErb et all[2010l : 
iFinkelstein et al.ll20lH )." In contrast the tremendous sen- 
sitivity of HST/WFC3 allows for > la detections of emis- 
sion lines to / > 2.5 x 10~^^ erg s~^ cm~^ for every ob- 
ject in a 2'.1 x 2'.1 field with 8 orbits. In addition, HST is 
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Fig. 1. — Line luminosity limits for detecting the Ho, H/3, [Olll], 
and [On] emission lines with 8 orbits of the HST/WFC3 G141 
grism (corresponding to a line flux limit of 2.5 X 10~^^ erg/s/cm'^). 
Dashed vertical lines indicate the redshift ranges at which each line 
is observed in the G141 wavelength range (1.1 < A < 1.7^m). In 
this paper we study galaxies with H/3 and [Olll] lines, selected at 
1.3 < 2 < 2.4. 

not limited by atmospheric emission or absorption in the 
NIR, and so spectral lines are detectable over the entire 
1.1-1.7 /im wavelength range of the WFC3 G141 grism. 
Figure [T] shows the detectability of emission lines with 
the G141 grism at various redshift ranges. 

Star formation processes were q uite different at z > 1 
than in the local universe (e.g. iPapovich et all 120051 : 
iReddv et al.ll2006l:lShim et al.ll2011[ ). and so it is not clear 
if the emission-line diagnostics used to characterize local 
galaxies can be directly applied to z > 1 systems. For 
example, the bulk oi z 2 galaxies have higher ratios 
of coUisionally excited to recombin ation lines ([O iii]/H/? 
and [Nil] /Ha) than local galaxies (jShaplev et al.ll2005"bl: 
lErb et al.l [20061: iKriek et al.ll2007D . Local galaxies with 
high SFR tend to have higher emission-line ratios than 
galaxies with lower SFR, suggesting that rapi d star for- 
mation leads to differen t ionizati on properties (|Liu et al.l 
[2008: Brinchmann. Pettini fc Chariot 200l . Since z>l 
galaxies typically ha ve higher SFR than lo cal analogs of 
the same mass (e.g. IPapovich et al.l [2005D . their higher 
line ratios ma y simply be t he res ult of their higher 
SFR. However I Wright et all (|2010D used spatially re- 
solved spectroscopy of one galaxy at z = 1.6 to show 
that these line ratios increase in the center of the galaxy, 
with normal star-forming (SF) line ratios at outer aper- 
tures. This suggests that weak or obscured active 
galactic nucleus (AGN) activity may be the cause of 
higher line ratios m. z 2 galaxies, perhaps because 
AGN activity is correlated with SFR as observed locally 
(jKauffmann et al.ll2003l:[Shi et al.,.20091 . There are many 
predictions for the density of obscured and Compton- 
thick (CT, Nh > 10^4 c m^) AGNs at L 5 < z < 3 from 
the X-ray background (jGilli et al.l 120071: iTreister et~all 
[200l and from IR -exces s galaxy counts and X-ray stack- 
ing (IDaddi et al l I2007t [Martinez-Sansigre et all 120071: 
iFiore et al.l 120091: IGeorgantopoulos et al l 120111 ) but to 
date only three C T AGNs have been ind i vidually con- 
firmed at z ~ 2 (lAlexander et al.l 120051: lErlund et"all 
[20081 1 Alexander et al.ll2008n . Obscured AGNs are missed 
by even the deepest X-ray surveys, but can still ex- 



hibit the narrow emi ssion-line signature of an AGN (e.g., 
LaMa ssa et al.|[20lll) . Measuring the density of obscured 
AGNs at z ~ 2 is especially important to reveal how 
much black hole growth occurs in an enshrouded phase. 

Here we use 8.5 orbits of WFC3 G141 spectroscopy to 
study the properties of emission-line galaxies at 1.3 < 
z < 2.4 in the H ubble Ultra Deep Field (HUDF, 
iBeckwith et"aIll2006D . These data were taken as part 
of the supernova follow-up program in the Cosmic As- 
sembly Near-Infrared Deep Extragalactic Legac y Sur- 
vey f CANDEL jg, Grog in et al . 2011; Kockcmocr'eFal] 
I201I . The supernova c andidate identified in the HUDF 
region is described by iRodnev et al.l (|201lD . and this 
work makes use of the spectra that were serendipitously 
obtained for hundreds of other targets in the 2.1' x 2.1' 
field. In particular we focus on the 28 galaxies exhibit- 
ing H/3 and [Oiii] emission lines in the observed grism 
wavel ength range (a t 1.3 < z < 2.4). Like the study by 
Wrig ht et al.l (|2010f ), the high spatial resolution of the 
HST/WFC3 similarly allows for studies of spatial gradi- 
ents in emission-line properties. Unlike an integral field 
unit spectrograph, however, the WFC3 grism has the 
unique ability to do this for hundreds of objects in the 
field of view simultaneously. 

In §2 we describe the observations, data reduction, and 
emission line luminosity limits. In §3 we describe our 
method for calculating emission lines, complicated by the 
limited spectral resolution of the WFC3 G141 grism. In 
§4 we stack the spatially resolved spectra to show that 
both the emission-line properties and the stacked X-ray 
data are suggestive of weak AGN activity in at least 
some of these galaxies. We summarize our results in 
§5. Throughout the paper we adopt a cosmology with 
h = 0.70, = 0.3, = 0.7. 

2. HST/WFC3 SLITLESS GRISM OBSERVATIONS 

The HUDF was observed with the HST/WFC3 GUI 
grism on 26-27 Oct 2010 and 1 Nov 2010 for a total of 8.5 
orbits. The data were acquired in two different orients, 
with 2.5 orbits at an orient of PA = 176° on the sky and 6 
orbits at an orient of PA = 182.4°. The data for ea ch ori- 
ent w ere reduced using the aXe software (Kiimmel et al.l 
120091 available at http : //ax e . stsci . edu/axe/[ ), which 
flux- and wavelength-calibrated the data and provided 
co-added two-dimensional (2D) spectra as well as opti- 
mally extracted one-dimensional (ID) spectra. The us- 
able wavelength range of each reduced spectrum is 1.1 < 
A < 1.7/im, with a resolution of i? ~ 130 (46.5A/pixel 
for a point source and O'.'13/pixel). 

The G141 grism is slitless, and this introduces a num- 
ber of complexities when interpreting the data. Spatially 
extended sources have features which appear extended 
in the wavelength direction, and so broad emission or 
absorption lines can be caused by either large velocity 
widths or by the size of the object. In addition, spec- 
tra frequently suffer from contamination by 0th order, 
1st order, and 2nd order spectra of nearby sources. For 
this reason we did not combine the two orients, instead 
analyzing each orient independently and using them to 
identify contamination. 

We calculated redshifts for all 455 targets detected 
in the HUDF G141 data. The best redshift was found 

http://candels.ucolick.org 
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using a cross-correlation template-fitting IDL algorithm 
based on the publicly available idlspec2d package writ- 
ten by D. SchlegeO Previous s pectro s copic redshifts 
JWolf et al. 2004; Mignoli_eLaD l2005t iVanzella et all 
[2(308; Straughn et al. 2008) and photometric redshifts 
IWolf et a l. 2004; Grazian e t.aL 200 6; Ryan et al.l 120071 : 
iWuvts et"a l. 2008; Cardam one et al.ll2bl0t ) were used as 
priors. We used two templates: a star-forming emission- 
line galaxy and a passive absorption line galaxy. Each 
template was constructed using the high-resolution op- 
tical (3800 < A < 9200A) templates from the SDSS li- 
brarjH combined with supplemental low-resolution UV 
and NIR data from the SWIRE photometric templates 
(jPolletta et al.|[2007[) , with the final templates resampled 
to match the spectral resolution of the WFC3 G141 data. 
All spectra and their redshift fits were visually inspected 
to ensure quality. 

Of the 455 targets, we selected the 28 galaxies at 1.3 < 
z < 2.4 and H/3 or [Om] emission-line luminosities of 
L > 1.8 X lO'*^ erg s~^. This luminosity limit represents 
the flux limit for detecting emission lines at the la level 
(/ > 2.5 X 10"^*^ erg s'^) at z = 2.4. For aU 28 galaxies 
we detect both H/3 and [Oiii] in both grism pointings, 
and as a result their redshifts are unambiguous. Table 
[1] shows the properties of these 28 emission- line galaxies. 
Their median redshift is 1.86. 

2.1. Stellar Masses and Rest-frame Colors 

Rest-frame magnitudes and stellar masses were cal- 
culated from the G OODS-MUSIC photometry catalog 
(jSantini et al.l 120091) by fitting a set of template spec- 
tral energy distributions drawn from the PEGASE stellar 
population models fsee iFioc fc Rocca-Volmerangelll997l 
for a description of an earlier version of this stellar pop- 
ulation model). We allow dust attenuation following 
ICalzettil (pOOL) . with values of E{B — V)gas between - 
0.05 and 1.5. (Slightly negative E{B — V) avoids biasing 
the mean properties of the population by accounting for 
the rare cases when photometric error gives fits that pre- 
fer slightly negative attenuation. In practice all of our 
emission-line galaxies have best-fit E{B — V) > any- 
way.) All templates have solar metallicity, and included 
a broad range of exponentially-decreasing, constant, or 
exponentially- rising star formation beginning at zj ^4. 
A grid of stellar population templates were compared 
with photometric data points of each galaxy, given the 
spectroscopic redshift from the grism data, and the best 
matching template was used to calculate rest-frame col- 
ors and a stellar M/L ratio. Compared to galaxies with 
independently-estimated stellar masses (from indepen- 
dent photometry of si milar but not ident ical datasets) 
in the GOODS-S field (IWuvts et al.l 120081) . our method 
yields similar masses for intensely star-forming galaxies 
both with and without dust attenuation (i.e., the galax- 
ies of interest here), with a scatter of 0.3 dex. Our masses 
are 0.1-0.25 dex larger for moderately star-forming galax- 
ies or non-star-forming galaxies because our templates 
include olde r stell ar populations than those used by 
IWuvts et al.l (|2008f ). Rest- frame colors are reproduced to 
within 0.1 mag in rest- frame {U — B) color for all galaxy 

http:/ /spectro. princeton.edu/idlspoc2d_install. html 
http;/ /www. sdss.org/dr5/algorithms/spectcmplates/ 
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Fig. 2. — Rest-frame {U — B) color and stellar mass for our 28 
emission-line galaxies (black crosses). Typical la errors are ~0.1 
magnitudes in {U — B) and ~0.3 dex in mass. For comparison 
filled gray circles show the 137 galaxies at 1.3 < 2 < 2.4 in the 
HUDF and gray contours and points show the entire set of 4140 
GOODS-S galaxies in the redshift range, both with masses and 
rest-frame quantities calculated from the same photometry and 
with the same method. Blue diamonds show a comparison sample 
of 2 ~ 2 galaxies from Shaolev ct al. (2005al): this is the parent 
sample of most previous studies of emission-line galaxies at z ~ 2. 
Arrows show the effects on {U — B) color if the H/3 and [Olll] 
emission lines are removed (although the [On] emission line, which 
is not measured here, could have the opposite effect on {U — B) 
color). The WFC3 grism reveals galaxies with masses a factor of 
~ 7 lower than previous studies of 2 ~ 2 emission-line galaxies. 

types. Full details on t he derived masses a nd rest-frame 



quantities are gi ven by iBell et ahl (lin prepD . 
As shown by Ivan der Wei et all ( 201 ID . strong emis- 



sion lines can have significant effects on broad-band col- 
ors in z ~ 2 galaxies. In our galaxies the median to- 
tal equivalent width (EW) of the H/3 and [Oiii] lines is 
EWabs = 220A: this brightens the rest-frame B magni- 
tude by 20% and reddens the {U - B) color by 0.2 
magnitudes. Despite these effects, we do not correct 
rest-frame magnitudes for emission lines for two reasons. 
First, the [On] emission line could have the opposite ef- 
fect, brightening U and causing U — B to appear bluer, 
but we do not measure this emission line because it is not 
generally present in our observed spect ral range. In ad- 
dition, the SF/AGN diagnostic we use (lYan et al.ll20Tll) 
was calibrated with [U — B) colors uncorrected for emis- 
sion lines, and so we do not correct our colors to remain 
consistent. 

Figure [5] shows the rest-frame {U — B) colors and 
stellar masses for our 28 emission-line galaxies. Also 
shown by gray filled circles are the 137 galaxies in the 
HUDF with photometric or spectroscopic redshifts at 
the same 1.3 < z < 2.4, with rest-frame colors and 
stellar masses computed in an identical fashion. The 
emission-line galaxies studied here are rather common at 
1.3 < z < 2.4, making up ~20% of the total H < 26 
galaxy population. Gray contours additionally show 
the set of all 4140 GOODS-S galaxies with photomet- 
ric or spectroscopic redshifts in the same 1.3 < z < 2.4 
redshift range, and our galaxies have a roughly simi- 
lar distribution of mass and slightly bluer colors than 
the larger parent distribution. For comparison with pre- 
vious sa mples, we also show the 50 galaxies at 1.3 < 
z < 2.4 of lShaple"v~et al.. (,2005ai) . using their masses and 



4 



Trump et al. 



TABLE 1 

HUDF WFC3/G141 Emission-Line Galaxies 
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(U—B) colors calculated from their ph otometry using the 
publi cly available kcor rect package (Blanton fc RoweisI 
I2007D . The galaxies of iShaplev et all (|2005aD form the 
parent samp le for most preyious z ^ 2 emission-line 
studies (e.g. lShapley et al.ll2005bl: lE7b et al.l[2006l[20Tol: 



IWright et al.i2010l ). Thel^edian mass of the grism galax- 
ies is M^, ~ lO^-^ Mfn, a factor o f ~ 7 lo wer than the me- 
dian mass of the IShapley et al.l (j2005aD sample. Several 
of the emission-line galaxies from the WFC3 grism data 
appear to be star-forming dwarf galaxies with masses of 
only - IO^Mq. 

2.2. X-ray Data 

The HUDF has the deepest X-ray data in the sky, 
with 4 JVIs of Chandra data. Imaging was obtained in 
54 observations (obsIDs) over the course of 3 Chan- 
dra observing cycles in 2000, 2007 and 2010 using the 
Advanc ed CCD Imaging S pectrometer imaging array 
(ACIS-l lGarmire et al.ll2003D . The data were reduced us- 
in g CI AO v4. 2 a ccord ing to the basic procedure described 
in ILaird et al.l (|2009D . Before combining the observa- 
tions, the astrometry of each obs ID was registered to t hat 
of the GOODS-IVIUSYC survey (|Gawiser et al.ll2006l) by 
matching the positions of bright X-Ray sources to H- 
band selected sources, using the tool reproject_aspect. 
Source detection was carried out according to the method 
described in [Laird et al.. (,2009i) . and none of our galaxies 



are detected. We instead assign Bcr upper limits calcu- 
lated using the background and PSF size from the loca- 
tion of each galaxy on the C/ianrfra/ACIS detector. 

3. LINE MEASUREMENTS 

Because the grism is slitless, any spatial extent along 
the dispersion direction will artificially broaden spec- 
tral features. The [Oiii]A4959 and [Olii]A5007 lines are 
often completely blended together in resolved galaxies, 
and larger systems additionally have H/3A4861 partially 
blended with the two [Oiii] lines. For this reason we 
do not attempt to individually measure each of the H/3 
[Oiil] lines directly from the data. Instead we subtract 
a linear continuum in the region 4770 < A < 5IOOA and 
fit the residual with a mixture of three Gaussians. The 
three Gaussians are constrained to be centered within 
20A (-1200 km s-i) of 4861 A, 4959A, and 5007A, in 
the rest frame. Because the two [Oiii] lines are blended, 
we measure the [Oiii]A5007 flux as 3/4 of the total [Qui] 
flux (i.e., 3/4 of the the sum of the two Gaussians cen- 
tered near 4959A and 5007A; IStorev fc Zeinro^ [200(11) . 
Errors in line flux (and /([Oili])//(H/3) and i[oiii]) are 
calculated by bootstrapping 10 000 realizations of the 
resampled data. 

Figure [3] shows the H/3 and [Qui] line regions for eight 
emission-line galaxies. The 2D dispersed image, ID spec- 
trum, and line fit is shown for each object. 
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Fig. 3. — Examples of fits to the H/3 and [Olll] emission lines for eight of the emission-line galaxies. Dashed vertical lines indicate the 
line centers for the H/3A4861, [Olll]A4959, and [Oni]A5007 emission lines. Because the G141 grism has such low spectral resolution, the 
H/3 and [Olll] linos are frequently blended. We subtract the continuum, shown by the red line, and fit the entire system with a mixture of 
three Gaussians, shown by the green line. The [Olu]A5007 flux is measured as 3/4 of the total [Olll] flux. 

The classical "BPT " diagnostic 

(jBaldwin. Phillips fc TerlevichI [TOSl uses emission-line 
ratios (namely, /([Olii])//(H/3) and /([Nii])//(Ha)) 
to distinguish galaxies dominated by star formation 
or black hole activity in the local universe (see also 
IVeilleux fc Ost erbrock 1987; Kewley et al. 2001). How- 
ever in the WFC3 grism data, Ha and [Nil] are 
hopelessly blended, and additionally are redshifted 
beyond the obse rved grism redshift range at z > 1.6. 
lYan et ahl (|2011| ) showed that it is possible to construct 
an alternate version of the BPT diagnostic which 
replaces /([Nii])//(Hq;) with rest-frame {U — B) color. 
This modified BPT diagnostic is especially useful for 
interpreting the WFC3 grism data. 



An alternative SF/AGN classification di agnostic possi- 
ble fo r our galaxies is the "MEx" diagram (jJuneau et al.l 
1201 If ), which replaces the {U — B) color with stellar mass 
M*. Stellar mass, however, is dependent on the physics 
of the model, and could be systematically inaccurate due 
to the constraints of available templates. Meanwhile rest- 
frame color does not depend on the physics of the model, 
only its goodness of fit to the data. F or this rea s on we 
prefer {U-B) over M*, noting that the lYan et all ()2011[ ) 
diagnostic agrees well with the "MEx" method (see Ap- 
pendix A of Jun eau et al.i,2011i) . 

Figure [3] shows our z ^ 2 e mission-lin e gala xies on 
the modified BPT diagnostic of lYan et al.l ()2011|) . along 
with comparison sam ples at z ~ an d z ^ 0.6 (als o 
from lYan et "aI|[20Tl and z ~ 1.2 f from ' Liu et "aI1l200l . 
Galaxies dwelling on the left and bottom of each panel 
tend to be SF dominated, while systems in the upper 
right are AGN dominated. Low metallicity galaxies tend 
to have bluer colors and higher /([Oiii])//(H/3) ratios, 
and indeed the locus of SF galaxies tends to shift to the 
upper left from z to z ^ 0.6 to z > 1. Our galaxies 



generally lie in same region as previously studied z > 1 
galaxies. 

Directly interpreting the position of z > 1 galaxies on 
the modified BPT diagnostic in Figure U is challenging 
for severa l reasons, however. First, the dividing line of 
lYan et al.l (2011) between SF galaxies and AGN was em- 
pirically calibrated at z < 1 only, and it is unclear if this 
dividing line changes with redshift. Local gal axies with 
higher SFR tend to have higher line ratio s ()Liu et al.l 
120081: IBrinchmann. Pettini fc CharlotI [20081 ) . and since 
z ~ 2 galaxies have higher SFR galaxies th an local galax- 
ies of similar mass (|Papovichet all [2001 . we might al- 
ready expect them to have a higher /([Oiii])//(H/3) ra- 
tios from both their rapid star formation and low metal- 
licity. On the other hand, active galaxies might also be 
more likely to be confuse d as SF-dominated in lo w-mass, 
low-metallicity galaxies. Ilzotov fc ThuanI (j2008D showed 
that AGN accreting as high as the Eddington limit con- 
tribute less than 10% to the emission-line ratios in metal- 
poor dwarf galaxies. Our z ~ 2 galaxies have similarly 
low masses, and the emission from their low-mass black 
holes would be overwhelmed by a SFR of a few solar 
masses per year. 

For these reason it is difficult to interpret the abso- 
lute position of objects, and more instructive to consider 
relative positions of the inner and outer regions of the 
galaxies, as well as differences in the profiles of [Oiii] 
and H/3. 

3.1. Stacked HST/WFC3 Data 

The slitless WFC3 grism provides spatially resolved 
spectroscopy, with a scale of O'.'IS per pixel. This is 
particularly useful for studying gradients of emission- 
line ratios: most local galaxies have a higher ratio of 
/([Oiii]) //(H/3) in outer regions due to negative metal- 
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Fig. 4. — Ratio of [Olll] to H/3 with rest-frame {U — B) color. The sohd line shows the division between AGN (upper right) and star 
formation (left and bottom) from Yan ct al. (2011). Crosses show the 28 2 ~ 2 galaxies from the WFC3 G141 data. Contours and points 
show comparison samples (also from Yan ct al. 2011) at z^' (left panel) and z ~ 0.6 (middle panel). In the right panel, gray diamonds 
show a comparison sample at z ~ 1.2 from Li u et aTl II2008I V The red points show the colors and line ratios of the stacked data: the 
diamond shows the inner 0'.'39 region and the circle shows the outer 0'.'39 — 0'.'78 region. The 2 ~ and z ~ 0.6 comparison samples have 
the same line luminosity constraint as the WFC3 G141 data, with either H/3 or [Olll] of L > 1.8 X 10*^ erg s~^. From z 0.1 to z ~ 0.6 
to 2 ~ 1.2, galaxies have lower metallicities and con sequently m ove to the upper left of the figure, and our 2 ~ 2 galaxies lie in very similar 
space to the previously studied 2 ~ 1.2 galaxies of ILiu et al.l | [200S ) . The line ratio and color of the inner region of the stacked data is 
marginally (I.Oct) more AGN-like than the outer region. 



licity gradients (e.g.. Ivan Zee et al1ll998| ). while a higher 
ratio in a galaxy nucleus would suggest the presence of 
an AGN. In general the signal-to- noise (S/N) of our data 
is too low to measure line ratio gradients for individual 
galaxies; instead we stack the 2D spectroscopy to study 
line ratio gradients. 

We stack the 2D images for each galaxy, equal- 
weighting each galaxy by normalizing its total (H/3 -|- 
[Oiii]A4959 -I- [Oiii]A5007) hue emission. Images are 
aligned by the peak of the emission (dominated by the 
[Oiii]A5007 line). This stacked spectrum allows us to 
study the spatial extent of the H/3 and [Oiii] emission 
lines. 

One-dimensional spatial profiles can be constructed by 
extracting each emission line along the cross-dispersion 
direction. Figure O shows the spatial profiles for the 
H/3A4861 and [Oiii]A5007 emission lines, along with the 
collapsed cross-dispersion profile of the entire 2D spec- 
trum. Each emission line is extracted from the cen- 
tral three pixels (140 A) about the line center. The 
H/3 line has a cross-dispersion full-width half-maximum 
(FWHM) of 7.2 pixels (0'.'94, ~8 kpc), while the [Oiii] 
hue is narrower with a FWHM of 4.9 pixels (0'.'64, 
~5 kpc). 

We test the significance of the [Oiii] line being more 
spatially compact than the H/3 line by bootstrapping 
100 000 resampled profiles. We construct each boot- 
strapped profile by adding a normally-distributed ran- 
dom error (with a given by the propagated error of the 
stacked profile) to the original data, and then measure 
FWHM([Oiii]) and FWHM(H/?). The fraction of boot- 
strapped profiles with FWHM([0iii])iFWHM(H;3) gives 
the significance of [Oiii] being more compact than H/3. 
This occurs in 98.1% of the resampled data sets, and so 




Pixel 

Fig. 5. — Spatial profiles of the H/J (blue line) and [Olll] (red 
line) emission lines, along with the entire spectrum (black line, 
offset by +0.4 in flux). Each ID profile is extracted along the 
cross-dispersion direction, using the three pixels (140A) about each 
line center for the emission lines or the entire wavelength range for 
the entire spectrum profile. The H/3 profile is more extended than 
the [Oiii] profile with 98.1% confidence, suggesting that at least 
some of the objects have a compact AGN contributing to the [Qui] 
emission. 

the [Oiii] line is more compact than the H/3 line with 
98.1% confidence. 

We also measure the spatially resolved /([Oiii])//(H/3) 
ratio and {U — B) in the inner and outer regions of the 
stacked spectrum. We designate the "inner region" as 
the central 3x3 pixel box (0'.'39 x 0'.'39, or 3.2 x 3.2 kpc 
&t z — 1.86) in the cross-dispersion direction, and the 
"outer region" as the two boxes 3 pixels wide from pixel 
numbers 4 through 6 (0'.'39 — 0'.'78) above and below the 
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Fig. 6. — X-ray luminosity with [Olll] luminosity for the 28 z ~ 2 
galaxies from the WFC3 G141 data. Since all of our galaxies are 
undetected in X-rays, we show them as X-ray upper limits. The red 
diamond shows the stacked X-ray data (and mean i[oni]) for the 
sample. The dashed line s hows the typical relat ion between Lx and 
^[Onil for Type 2 AGNs ^LaMassa et al.H2011D and the dotted line 
shows the typical relation for star-forming galaxies IRovilos et al.l 
120091 ). both at 2 ~ 0, with the Icr scatter about each relation 
indicated by the gray shaded area. The position of the stacked 
X-ray data in hx — i[oni] space suggests that an AGN is present 
in at least some of the galaxies. 

center in the cross-dispersion direction. Line ratios are 
calculated for both the inner and outer regions, using 
the same method of continuum subtraction and Gaussian 
fitting described above. 

Color gradients were calculat ed using the CANDEL S 
WFC3 y, J, and H imaging (jKoekemoer et al.ll20Tll ). 
The WFC3 images are drizzled to a finer resolution 
(0'.'06) than the WFC3 spectroscopy, and so we adopt 
0'.'36 (6 pixel circular diameter) as the inner region and 
0'.'36 — 0'.'72 (circular annulus from 6 to f 2 pixels) as the 
outer region. Rather than stacking, we measure aperture 
y, J, and H for each individual galaxy. We then treat 
the median Y — J (for low redshift 2; < 1.8) or J — iJ 
(for z > 1.8) difference from the inner (< 0'.'36) to outer 
(0'.'36 — 0'.'72) aperture as the change in [U — i?)iost over 
the same aperture. 

As Figure S] shows, the stacked inner region has 
marginally redder colors and slightly higher line ratios 
than the outer region. In other words, the inner region 
is more like an AGN than the outer region. However the 
large errors (particularly on the outer ratio) mean that 
this is significant only at the I.Oct level (67%). This is 
probably because the 0'.'39 size of the inner region corre- 
sponds to the inner 3.2 kpc of the galaxy (at the median 
redshift z = 1.86), and so includes much more than the 
typical ■^l kpc narrow line region (NLR) of an AGN. The 
marginal difference between inner and outer line ratios is 
merely suggestive, but is reinforced by the far more sig- 
nificant (98.1% confidence) result that the [Oiii] profile 
is more compact than H/3. 

3.2. Stacked Chandra Data 

In Figure H] we show the full-band (observed 0.5- 
10 keV) X-ray luminosity upper limits with the 
[Olil]A5007 emission line luminosity. The dashed line 
shows the typical Lx — ^[oiii] rela tion for [Oiii]-selected 
narrow-line AGN ([LaMassa et al.|[2011|) , with associated 



Ict scatter of ^^0.5 dex indicated by the gray shaded area, 
and the dotted line shows the typical Lx — i]oiii] rela- 
tion for star- forming galaxies ()Rovilos et al.|[2009l) . with 
the gray shaded area similarly indicating the Ict scatter 
of ^0.3 dex. 

We stack the X-ray data for the undetected sources to 
learn whether the population has Lx — ^[oiii] dominated 
by black hole activity or star formation. The stacked 
data is marginally detected (2.9ct significant) with a flux 
of /0.5-iokoV = 1-1 X 10~^^ erg cm~^. At the me- 
dian redshift of our sample, z = 1.86, this translates to 
an X-ray luminosity of io.5~iokcV — 2.7 x lO'*^ erg s~^. 
Because the stack is only marginally detected in the full 
band, we cannot comment on its hardness. Unless the 
typical Lx — ^[oiii] relation for star formation changes 
with redshift, the position of the stacked data is interme- 
diate between the loci of AGN and star formation. This 
suggests that at least some of the galaxies have an AGN 
contributing to the [Olll] and X-ray emission. 

4. NATURE OF Z - 2 EMISSION-LINE GALAXIES 

The stacked set of our z ~ 2 galaxies shows a more 
compact [Oiii] profile than the H/3 emission line. This is 
the reverse of w hat is expected for ty pical galaxy metal- 
licity gradients (jvan Zee et al ] [l99l . but is expected if 
an AGN in the galaxy's center contributes to the line 
emission. We find a marginally higher [Qui] /H/3 in 
the center of the st acked data, s i milar to the higher- 
confidence result of iWright et al.l (|2010( ) , who found a 
steeper [Oiii]/H/3 ratio in the center of one z = 1.6 
galaxy. In addition the stacked Chandra data suggests 
an AGN component in at least some objects, since the 
position of the stacked data in Lx — ^[oiii] is inter- 
mediate between the lines of star formation and black 
hole activity. For these reasons we suggest that some of 
these z ~ 2 emission-line galaxies, and tho se with similar 
line ratios and colors from other authors (Shaple v et al.l 
iJOOSb; Erb et al. 2006; Krick ct al. 2007), probably con- 
tain weak AGN in their centers. 

Because we rely on stacked data to reveal an AGN, it is 
difficult to know the frequency of nuclear activity in these 
galaxies. In the local uni verse, AGN act ivity is known to 
correlate with SFR (e.g. iShi et alll2009| ). so AGN activ- 
ity might be expected in our galaxies given their bright 
emission lines and presumably high star formation rates. 
While some authors have argued that AGN are rare in lo - 
cal galaxies with Af* < lO^^M© (jKauffmann et al.ll2003f l. 
recent work suggests instead tha t AGN incidence is con- 
stant with host mass at z < 1 ()Aird et al.ll2011h . Our 
sample has a median mass of = 10^ -^Mq, and so 
unless a few massive objects are dominating the stacked 
X-ray and grism spectrum, this suggests that AGN are 
common in low-mass star- forming galaxies at z ^ 2. 
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